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By Karl Marguerre 



The present report is a very simplifying derivation . 
for the results of an investigation entitled: "The Load 
Capacity of a Plate Strip Stressed in Compression "beyond 
the Buckling Limit" and published in 1937 in the Zeit- 
sohrift fur angewandte Mathematik und Mechanik. 

Following the discussion of the methods and results 
of other authors, the writer suggests an extension which 
is very desirable from the point of view of airplane de- 
sign problems. It affords a practical theoretically 
evolved formula for the apparent width under an apprecia- 
bly exceeded buckling load. 



SUMMARY 



The present extension of the customary stability in-' 
vestigation to include the supercritical range, proceeds 
in two steps. The first step considers the buckling form 

TT3C TTTT 

w = f cos -=— cos -— known from elementary theory, pre- 
l> b 

serves the higher terms in f and yields, with the aid of 
the principle of. virtual displacements, a relation \?hich 
gives the decrease of the apparent strain stiffness at the 
instant of buckling (analytically expressed "the tangent 
to the now stress-strain curve above the critical load" 
(equation 5.5)). 

The second step evolves on the basis of a formula 
containing several arbitrary values, from which the prob- 
ably produced buckling form; with a greatly exceeded crit- 
ical point, can be computed, and which affords a stress- 
strain curve (fig. 2) which reproduces with sufficient 
agreement the actual Conditions existing in a zone 



*"Die mittragende Breite der gedruckten Platte." Luft- 
f ahrtf orschung, vol. 14, no. 3, March 20, 1937, pp. 
121-128. 
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c* < e x < 20"€* (as proved "by comparison with experimen- 
tal results). The apparent width is conveniently calcu- 
lable (fig. 3) within the cited range with the aid of for- 
mula (7.8) cr (7.9), which approximately comprise the re- 
sult of the theory. 

I. INTRODUCTION 



The following investigation treats the load capacity 
of a rectangular plate stressed in compression in one di- 
rection (x) beyond the buckling limit. The plate is ro- 
tatably (i.e., free from moments) supported at t.11 four 
sides by bending-resi stant "beams. 

Before buckling, the axial compression p x = - a x is 

uniformly distributed and proportionality exi sts .between 
the crushing €i = - £ and the compression p x (accord- 
ing to the law of elasticity). Above the critical value 
€*, of the crushing, the strip buckles - more in the mid- 
dle than near the restrained sides - resulting in nonuni- 
formly distributed axial compression; the centroidal axes, 
as commonly expressed, "do no longer fully contribute." 
The sought-for factor is the condition of form change in 
the buckled sheet and in particular, the new stress-strain 
curve; i.e., the relation between the mean value: 

?X " b / »X ^ = P ! 
A/2 

of the compression and the mean crushing Ci (the crush- 
ing of the longitudinals): 

P x = P x (€ x ) for t x > e* 

The crushing e x (respectively, the amount of the 
pushing together of the transverse beams e i I ) is chosen 
as the first independent variable of our problem; whilo as 
second independent variable, the lateral displacement 
(pushing together) e 2 h of the longitudinal sides or 
else the mean value of the compression in transverse di- 
rection is introduced. 
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II. THE FUNDAMENTAL EQUATIONS 

The hypotheses of classical plate theory (preserva- 
tion of the normals, etc.) allow us to express the total 
stress-form change condition of the thin plate as func- 
tion of the three displacements u = u (x,y), v = v (x,y), 
w = w (x,y) of the plate middle. The problem of plate 
stretch (displacements u,v) may be reduced to the bipoten- 
tial equation A A $ = 0 by having recourse to a "stress 
function" <I> ; the deflection w follows the "plate equa- 
tion 11 - i.e., the bipotential equation with interference 
term. The equations for $ and w are unrelated* 

Premise of this theory is the fundamental assumption 
of the "linearized" elasticity theory: that all displace- 
ments relative to the dimensions of the body ,, particularly 
as regards the plate thickness s, are small. 

But a thin-walled sheet may undergo elastic deflec- 
tions amounting to multiples of its plate thickness; for 
the treatment of problems of that kind the linearized 
plate theory falls short. Now, a very practical and em- 
pirically closely agreeing theory is arrived at "by extend- 
ing the elementary formulas so that the quadratic portions • 

in the deflections w are retained in the changes of the 
coefficients of the linear element ("strain" and "slip-> 
page") while, as before, the higher powers of the stretch 
u,v and likewise the products s w and z w (because s 
is of the same order of magnitude as w) are stricken from 
higher than the second order. 

If Y denotes_the changes of the coefficients of the 
line element, and Y its mean values over the plate thick- 
ness (strain portions), we have: 



A more elaborate argumentation of this hypothesis is 
found in a report by E. Trefftz, entitled "On the Deriv- 
atives of Stability Criteria ... III." Intern. Mech. 
Kongress, Stockholm, vol. Ill, 1930, p. 44; s.a. Handbuch 
der Physik, vol. VI, p. 56. 
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= Y" ia - 2 2 



3 s w 
3xdy ' 



2 £ + (1?) ■ 



'is 



du .[_ 8v + 8w 8w 
8y 8x 8x 8y 



(2.1) 



whereby the terms w x 2 t 7^2 j Wx Wy are new compared to 
classical theory. 

All other assumptions of th.e elasticity theory of 
small deformations may be retained unaltered. In particu- 
lar, Hooke's law in original form: 



%-»•(=¥■ + » ! ! 3 -). «v - 



E 



and the terms: 



- 1 [„ Yl1 



+ a. 



23 



y 2 



— + T Y 
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= ~ [(\i +Y 32 ) 2 - 2 (l-i>) (Y^Ym-Y^ 8 )]. j> 
1 



8 

23 [(cx x + ^y) 8 - 2 (l+v) (a x a jr -T 2 )] 



(2.3) 



for the form change energy (FB, abbreviated) a z P er unit 
volume retain their validity. 

The EE per unit surface a 2 is obtained, for exam- 
ple, from equation (2.3) by integrating z from - s/ 2 
to + s/2 as the sum of the "strain energy" a s and the 

"bending energy" a 3 

f 

a 2 = a 2 + a 2 



with 
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- 2-U-W) [(u x + | w x a ) (v y ■+ | w y a) 
- (u y + v x + w x w y) S J} 



For the practical calculation it is of advantage to 
introduce the stresses rather than the form changes in the 
strain portion a 2 . Dividing the stress a (x,y,z) pre- 
vailing at any point within the plate in the conventional 
manner into strain a and "bending stresses q (say, in 
form of 

Ox = 5* - Si , 0-y = CTy - ffy , T = T - 2* T ) 

the strain portion expresses itself through the stress 
mean values <r in the form of 



= p! + ^ - 2 (1 - v) (a x <J y - r ff )] (2.4 3 ) 



On the other hand, the strain stresses follow Cauchy's 
equilibrium equations: 

dcF x 3t~ 3t 3a v 

and the form of these equations allows us to replace the 
three unknown func 

tion" $ through: 

^ = 3>yy, Oy = T = - $ x y . (2.S) 

Hence the expression (2.4 a ) for the strain portion of FE 
"becomes : 

*a = 2 S e t(^x^yy) 2 ^ U+ ( *xx*yy-*X7 3 (2-* S ) 



three unknown functions a_, 0" , t "by one "stress func- 

x y 



and the total form change energy A stored in a plate of 
length I and width "b may, according to (2.4 3 ) and 
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(2.^) "be written as: 
1/2 fc/2 

-1/2 -.-b/2 



($ $ -$ 2 )"| 
v *xx*yy *xy ' J 



JLLjSi 
24 



C(Aw) S -2 (1-v) (w xx w y3r -VT xy 2 ) ] jdx dy 



(2.7) 



An equation "between stress function $ and deflection 
w is obtained if Hooke's law (2.2) is specially rewritten 
for the mean values (strain portions) of the stresses and 
strains: 



= E ' (u x + v v y + I (w x a + U Wy 2 )^) 
= E ' ( v y + u ^x + * (V + u w x 2 )) 



$ xy = G + u y + w x w y) 



(2.8) 



and the displacement u, v is eliminated from these three 
equations (reference l): 



A A $ = E (W x y 2 - W XX Wyy) 



(2.9) 



With the aid of equations (2.6) to (2.9) the condi- 
tion of stress and strain in the plate can be progress- 
ively determined (Ritz's method). 2 



III. DETERMINATION OP STRESS FUNCTION $ AND OF 
DEFLECTION f AS FUNCTION OF e a AND Cj, 



In contradistinction to the pure energy method em- 
ployed in the report quoted at the beginning, which makes 
explicit use of only equations (2.2) and (2.4i) we apply 



There is little hope for an exact integration of (2.9) to- 
gether with the "extended, plate theory" which might "be 
added as second equation to determine $ and w (cf. K. 
Marguerre: Z.f.a.M.M. , vol. 16, 1936, p.. 353), "because 
both equations are nott-linear in addition to being coupled 
(in contrast to the elementary theory). 
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a, "mixed" method. The deflection w is again expressed 
"by a Ritz formula (containing arbitrary values), hut we 
first compute. 3> from the differential, equation (2.9), 
and then only proceed to the energy expression in the form 
of (2.7) in order to determine, "by virtue of the minimal 
requirement of the principle of virtual displacement, the 
free values as functions of the independent £i and e 2 . 

Even though this method does not suggest itself as 
readily as the previously employed one, it has (and this 
is of particular advantage in more complicated cases - re- 
strcint, shear) the advantage of minimizing the paper 
work. 



w : 



In fact, the simplest Ritz formula 4 imaginahle for 

ttx Try 

w = f cos — cos ~ b~ ) (3.1) 

tfThat is, assuming that this form of "buckle produced at the 
instant of "buckling, is preserved in form even "beyond some 
distance after exceeding the "buckling load (i.e., that on- 
ly the free value f changes as the load increases), gives, 
for the right-hand side of equation (2.9): 



2tt y\ 

cos ! 

b / 



, TT f / 2TT x\ / 

- (l + cos Sp) (l + cos ^Z)] 



and a particular integral of the equation: 

tt 4 f 2 / 2tt x 2tt y 

g- o" (COS + COS 

2b I n I T> 



TT I / 2TT X £-TT y \ 

A A $ = - E — 3- — a ( °os + cos ) (3.2) 

2b I N 1 b / 



l Both methods are identical in nature "because the differ- 
ential equations (48) for the displacements u, v are sim- 
ply the equilibrium equations (2.5) written in the displace- 
ments. 

If I ^>b this theorem must he extended to include the 
nodal points in transverse direction. Then I denotes 
the distance of two nodal lines. 
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can "be given at once in the form of 

„ f 2 f/l\ 3 2tt x /b\ S 2tt y 1 

$ p a - E 3i Lw cos — + Cr ; cos — J 

The complete solution for $ satisfying the limit- 
ing conditions is obtained with appropriate solutions of 
the correlated homogeneous equation: 

A A * = 0 (3.3) 

Assuming predetermined displacements, the required 
limiting conditions are: 

u (± |, y) = T c x 1, t («. ± |) - J «. | (3.4!) 

for the displacements normal to the sides; 

▼ (± |, y) = ? e 2 y, u (x, ± |) = + £ t x (3.4 a ) 

for the displacements tangential to the sides, 

These conditions are not exactly satisfied "by combina- 
tion of the particular integral with a. finite number of so- 
lutions of equation (3.3). 

On the other hand, it has "been proved that the careful 
compliance of the limiting conditions for the displacement 
tangential to the sides has no marked effect on the sought- 
for stress and form change condition of the plate and par- 
ticularly, on the stress-strain curve. The requirement of 
vanishing derivation of normals instead of (3.4a) 



3v (± |. y) 3u (x. ± |) 



Si - »• —gr— - ° (3-4=) 

along the cited sides - mechanically expressed "disappear- 
ing shear" 5 - may in consequence be introduced. 



G \ 2' y / Lax dy 9x 6y 

1 T ( Xf ± = T|v + |u + |w |w 
G \ Zj j_9x 3y ax By, 



= ±1/2 
^r=±b/2 



Both expressions disappear, in fact. 
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The limiting conditions (3.4*) and (3,4 3 ) c an 
"by a happy chance, he complied with in very simple manner. 
Putting 



- $ = E 



fG) ! 



32|_ \ 



cos 



2tt x , fb\ 



2rr_y 1 
I J 



that is, 



^y si 



+ i r» x 2 + ^ p y ; 
2 2 2 1 



ni-|i C os ^-s . p 

8h 1 



- P x - *yy = ■ s g-f" c °* — - p * 



equation (2.8), that is: 

E u = $ - v 3> - 5 w 3 , E v_. = 

x yy xx 2 x 7 x 



(3.5) 



(3.6) 



" u $ yy - I V 



and (3,4 1 ), that is: 

n l ! 2 



EC! = - ^ 



($ yy - v $ xx ) dx - 



E 



,1/2 



w x a dx 



] 



a .pa 



- E 



TT a f 
8? 



2 Try „ rr 2 f 2 

cos — + p x - u p s + E --a— 



C 



2Try\ 

1 + cos -— J 



,V2 



E e 2 = - 



» I, 



2 h/2 
C*xx - v *yy> ^ - I ] V dy J 



8h a I 

give for the two constants p x and p 2 the equations: 



Pa = E» [e E + v € a 



(3.7) 
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while the condition - $» v . = T = 0 is identically ful— 

filled in the entire plate and consequently, also, in th e 
sides. 

The mechanical significance of constants p x and p a 

appears from the equations (3.6). Integration over the 
width and length of the plate gives: 



1/2 






T)/2 






I 


$ xx dx 


= ~ Pa ^ ' 


y $yy 


dy = - 


Pi * 


-1/2 

or 




h/2 


-h/2 


1/2 






« - i 


/ Px dx - 


-i 


A' *y 


dy 






:{/ 2 









In this manner the equations (3.7) give at once the sought- 
for relations existing between the mean values 'of the 
stress and the mean crushing €i and e 2 , leaving the do- 
flection f as the sole unknown quantity. 

The determination of f follows according to the 
principle of virtual displacement from the premise that 
the potential energy (i.e., the difference between the EE 
and the potential of the external forces) becomes a mini- 
mum. Having assumed rigid side beams while investigating 
the equilibrium condition by predetermined displacements 
of the sides, the external forces contribute no work on 
the sides (fixed during the virtual displacement), and the 
minimum requirement for the potential energy reduces to 
the minimum requirement for the EE itself. 

As a result, wo need to formulate only (cf. equations 
(2.7), (3.1), and (3.5)): 

A = / E (A$>2 - 2 (1+u) ($ *x ^y-V^ dx dy 

+ "l!" y'C(Aw) 3 - 2 C 1 -^ ( w xx w yy- w xy 2 )3 dx d ^ 

- sbl Fe (ulll + rr 4 f 4N \ . 1 , s a 2 (1+v) „ 1 

+ E« s 3 bl \ulf + rr 3 f ] _ a ?\ 
+ ^ + .^j _ A (p,f ) 
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\ 



and to substitute the strains for the stresses "by me*ans of 
43.X),,. giving:. ^ 

A=A(c,f)E»shl {-f + + V cxe 8 - f^f- [c s + i>) 

■t.«.(i*^)]*^[(^(^*i^}S]- 



and £^ii-»£l =' 0 gives for f the equation: 
3 f 



TT 



2 o 2 



12 b 



uL.il 

16 b 2 



[(3 - u 2 ) + l) + 4u $1 



(3.9) 



The three equations of (3.7) and (3.9) represent the sought* 
for stress-strain law above the critical load. 

The result is discussed in greater detail when con- 
fined to the case of the square plate (1 = b ) . 



IV. STRESS-STRAIN EQUATIONS FOR THE SQUARE PLATE 



The strain stresses are, according to (3.6) and (3.7) 



y 



= Et ( £l + V e 8 - n£l| (i + V ) 
<■ 8 b 3 

[l + (1 - v) cos 
= E* |e 2 + v € 1 - n£f£ (i + v ) 

L 8 b 3 



[ 



1 + (1 - u) COS 



2ttx 



} 



T =0 



(4.1) 



(4.2) 
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The shear stress disappears, the normal stresses are dis- 
tributed at right angles to the direction of action, ac- 
cording to a cosine law. 



For the mean values it gives with (3.9): 

n£f£ . f (e + € j . nisi] 

8b s 3 - V L 1 + v 3b s J 

»» ■ »■ Hi {««• - (l - B > £ = + rhj »iSi} 



3b 2 / 



(4.3) 



(4.4) 



For the maximum value of the compression p^ (occurring on 
the sides), it is: 

»b - «• Hi { <3+b) + v * HH]} (4 - 5) 

while for the minimum value (in the middle), we obtain: 



1 - u 



3 - v v. 



I (1-v) 



Cl 2£a + 1 - » 3b* 



) (4.6) 



In many cases it is of advantage to treat the stress 
and strain condition as function of the mean compression 
p x and p 3 ■ exerted on the side beams (as independent vari- 
able). The solution of the system (4.4) gives: 



S £j = 2pi + (1 - ») P 2 - P 0 | 
S «a = 2p 3 + (1 - V) Pl - Pq J 
where p„ is an abbreviation for 

* o 3b 
then (4.3) is replaced by 



(4.7) 



(4.8) 
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and (4.1) reduces to the simple form of 



P x = Pi - (Pj. + P a * - P Q ) cos 



2-rry 



Py = P 8 - (P X + P S " P 0 ) ^ 



(4.9) 



V. TEE SPECIAL CASES p = 0 AND _e 2 = 0 



The equations (4.1) to (4.9) "become considerably more 
simple when further divided, for two very important prac-' 
tical cases: 



ee -T-y sh if t i ng 
1 9 ng itudinal sides 
Po = 0 " 



M on shifting 
longitu di nal s ides 



From the second equation (4.4) and (4.7) follows: 



3 a 3 



TT e S 



e = l=u € 1_ 

2 2 1 2(l-u) 3d 3 " 

and consequently, 



i=2 n + i Po 
2 - i 2 ° 



8d^ 



1 / 1 TT 2 S 2 \ 

2 V 1 " 1-u 3 ~3-b 5 /' 



TT3f2 1+U / P Q \ 

(5.1l) 

from (4.3) and (4.8). 

On the other hand , (4.8) and (4.3) give at once: 

•ri TT 2 f 2 „ , \ TT 2 f 2 _ 1 (c » 1 TT 2 S 2 ^ 

(5.1 S ) 

With f = 0 these relations give the critical values: 



El ( f=0 ) = e* = —^-5 

1-v 3h 



Pl (f=0) = p' 



__1 nfsj = E e * 

1-u 2 3t> 2 



£l (f=0)= ... = _i_ 
Pi (f=o)= P ** = \l- 



_ e« -— nisi: _ E t e ** 



l+V 3d 3 



(5.2) 
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in accord with known results. 

The second equation of (4.4) discloses that the A pr'e"- 

fix for the transverse strain (or transverse compression^) 
changes when 

.a 




n - 3 b 



Referred to the critical crushing €* or e**, it means 
that (for a transverse contraction figure u = l/3): 

The transverse elongation The transverse compression 

"becomes a contraction becomes a tensile stress 

"by twice exceeded "buck- "by three times exceeded 

ling crushing buckling crushing 

With the abbreviations (5.2), the equations (4. 3) and (4.8) 
read: 



TT 2 f 2 



4V* 

E H ~J = 2 (Px - P*) 

4b 



4b 3 - u 1 



(5.4) 



Eliminating f from these two pairs of equations leaves 
the stress-strain law in the form: 



Pi 



- p* = 1 E (€!-€*) 



2 



Px - P 



(5.5) 

Note that the relationship of stress and strain 
remains linear even above the buckling load within the 
scope of our approximation; the "atmarent strain stiffness" 
is 

Unrelated to v for v = i 

3 

is exactly half as great as in the proper elastic range be- 
low the critical load. 

Observing that in both cases the compressive stress 



P L = E € X 
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prevails in the longitudinals, equation (5.5) gives for 
the apparent width "b m of the sheet, according to the 



equation 



- D m = 



IT 



\ 0 - ? --) 



\ (l + ii) 



b 



( 



1 + 



(3-v) 



lii^l p** 
~2" p L 



16 V 9 Pl y 



(5.6) 



The transverse tension set up by nonshifting sides 
causes a slight increase in apparent width (i.e., in the 
load capacity of the sheet) relative to the case of disap- 
pearing mean transverse stress. The same result obtains 
from a comparison of the minimum compression values pro- 
duced in the centroidal axis: Substituting the side 
stress p R and the critical compression p* and p** 

for e 1 and e* (e**), according to (4.5) and (5.2), 

equation (4.6) gives for p M the values: 



PM 



p M 



+ 2_Xl±2i v ** 



1 - V 



) 



= 1/5 p R + 4/5 p** (u=l/3) 

(5.7) 

where it will be observed that in the first case the cen- 
troidal axis rejects every compression rise above buck- 
ling, while in the second case, because of its better sup- 
port by the transverse fibers, it takes up (even though 
small) a part of the compression. 

The result expressed by equation (5.5) - that the ap- 
parent strain stiffness in a square plate at the. instant 
of buckling reduces by half, is quite generally valid for 



p 3 = const 4 0 



e 2 



= const 4 0 



The critical valtie itself is affected by the magnitude of 
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p 2 or € a $ 0. We find: 



'kr 



= p* - 



B € kr = P* - (1+u) P a 

= P kr - v (p* - p kr ) 

and from (4.7) and (4.4) follows: 
E (€i - € kr ) 

= [2 Pl - p kr - u(p*-p kr )] 

** [P kr - u(p*-P kr )3 



c kr = e** - c a 
Pkr = Si [€** - (1-v) c 3 ] 

= E» [e kr + u (e**-£ kr )] 



(5.8) 



= 2 <Pi - Pkr> 



Px - Pkr = E * 



1 - u 
3 - v> 



{[2.x ♦ (1 - U )c kr+ f^ ...] 
. [(3-p) £kr - JLggBi f] 



that is, precisely the law (5.5) except for the values 
Pkr"*"* according to (5.8) instead of the values p*, 

etc. 

The linear aspect of the stress-strain curve above 
the critical load is, of course, a result of the limita- 
tion to the approximate formula (3.1). As a matter of 
fact, equation (3.1) is valid only "at the very first in- 
stant" after exceeding the buckling load. The straight 
lines (5.5) are the tangents to the actual stress-strain 
curve, which, starting from these tangents, deflects down- 
ward (fig. 2). 



VI. COMPARISON WITH OTHER TEST DATA 



The investigations discussed thus far differ from 
older reports on this subject (Schnadel, referenced; Cox, 
reference 3, Yamamoto and Kondo , reference 4, Timoshenko, 
reference 5), in sofar as, other than hypothesis (3.1) 

TTX TTV 

w =• f cos — cos — (square plate) (6.1) 
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concerning the "buckling form, no other arbitrary assump- 
tions are introduced. 

Instead of determining the horizontal displacements ■• 
u and v explicitly - or by introducing the stress func- 
tion $ directly - the first three authors cited, intro- 

duce the assumption T = 0 or in other words, ^— = 0, 



6ay 

ten in the form: <7 



dx 



sr - °- 



By virtue of this assumption, 

X = j f ^ 

Hooke's law (2.8): '"^i ' 2 



6 a. 



may he writ- 
dx and the extended 



b/2 



b/2 



1 

b 



-b/2 

E 



dx = 



b/2 



E» 

2b 



-b/2 



(Vn + V V a2 ) dx 

b/2 



"^b/2 - y b/2 
changes in the special case of u = 0 to 



Cfc = J C(u(b/2) - u(-b/2)]+ E 



= E 



,V2 



- £ 



2b / x 
-b/2 

and correspondingly: 

,V2 



2b 



J 



w x 3 dx 



-'b/2 



(6.3) 



E 



c 2 + 



2b 



w. 



dy 



Admittedly, Schnadel adduces this relation by the stress- 
function method but the inference is right only for each 

single term of the formula w = 1^ f^, cos -Sjp cos ■ 

As soon as several summands are admitted simultaneously, 
the la.w of superposition naturally ceases to hold and, in 
fact, the mixed terms ^±^y z lead to a nondi sappearing 

shear stress whose omission, to be sure, renders the cal- 
culation much easier. 



is 
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Then the strain portion It of FE "becomes in the case of 
disappearing transverse contraction: 



I = 



s 

2E 



Es 
2 



JJ (^ 3 +a y 2 ) dx dy = U[ y^a x s dy + J* a/ dx] = 

) dy + f\e s /\ y8 dy) dx] 

'J 'J (a , 



(6.4) 

which means that this portion (by given Ritz formula 
for w) can he determined in very simple fashion. 

The energy of bending is given unchanged according to 
(2.4) or (2.7). 

Cox and Yamamoto-Kondo , after him, carry this sim-r 
plification one step farther by writing a (x,y) = 0 for 

the whole plate; that is, neglect even the work of strain 
performed by the transverse stresses. Schnadel, on the 
contrary, uses (6.4) as basis . Allowance for the trans- 
verse contraction is in none of the formulas possible. 

Cox's assumption Oy (x,y) = 0 contains two state- 
ments: that the transverse stress within is invariable 
/3a v (xy) \ 

( — £ = 0) and that it disappears at the side 

V By / 

(dy (x, ± b/2) = 0). In the practical, most important 

case of strong - i.e., especially undef 1 ectable longitudi- 
nals - the limiting condition is certainly not complied 
with; and so for this reason alone a direct comparison of 
Cox's results with ours is impossible^ Granted even that 
we visualize the limiting condition Oy = T = 0 to have 

been realized by some appropriate test arrangement, the 
first statement regarding the strain condition within, 
still remains fundamentally inadmissible. Because, in or- 
der to apply Ritz's method in a mathematically unobjec- 
tionable manner (that is, for example, preserve the known 
dictum that the true load capacity must lie below that 
computed by approximate formula), it is not permitted to 

7 

Since Schnadel proceeds from a somewhat different concept 
of the buckling process (division in internal and exter- 
nal energy), the true facts of the case are not apparent 
at once. 
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mix geometrical statements (concerning w) with mechanical 
ones (about t). Unless the mechanical theorem T = 0 
happens to "bo fulfilled by itself, the ■ ■ FE-. must ho re- 
duced by this assumption: concerning the direction of tho 
discrepancy from tho true valuo (let alone an estimation 
of the errors), nothing more can be said. This drawback 
is felt most when several approximate solutions are to be 
compared; because then it is impossible to deduce that the 
formula giving the lowest load capacity must be the better 
one. 

Mathematically, the omission of a part of the FE 
fares no better. So, for instance, Oox obtains for the 
stress-strain law above the buckling limit the expression 

Pl - p* = | (£i - e*) (6.5) 

i.e., a drop in (apparent) strain stiffness to one-third. 
But, if '«e compute the case of vanishing side stresses 
cfjr and t "exact" conformable to the method cited at the 
beginning of this report, we find: 

p x - p* = k E (£i - e*) (6.6) 

where the coefficient k ranges between 0.41 (u = 0) 
and 0,34 (u = •§ ) depending on the amount of transverse 
contraction. For v = 0.3 xz has the value 0.38; that is, 
below 0.50 (see equation (5.5)), but still noticeably above 
0.33. Consequently, the omission of the share of the 
FE originating from the transverse and the shearing stress 
is not without some effect. 

Cox's experimental results are in contradiction with 
his formula. (6.5), but confirm our formula (5.5) very sat- 
isfactorily. 

From Schnadel < s data the results (for p g =.0) given 
in the preceding sections can be deduced in complete agree- 
ment, because the shear stress actually does disappear within 
the validity range of w = f cos ~ cos ~, and the effect 

■ ■ 

of the transverse contraction cancels out in the case of 
displaceable longitudinal sides. 

■ Timosheiiko elects to proceed from the exactly valid 
expression (2.4j) rather than (6.4). But his calculation 
differs from the one given here by the introduction of 
approximate assumptions, each containing a free valuo for 



I 
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the displacements u and v (independent of the tt for- 
mula). The results reveal the strange fact that the stress 
maximum docs not occur on the side hut on the inside at a 
distance depending on the degree of loading, so that the 
stresses actually decrease even toward the side. Apart 
from that, there are (as in Yamamoto's case) ten sile stress- 
es in the center of the hay under greatly exceeded "buckling 
load. Instead of the formulas (5.5), through which the de- 
crease in apparent strain stiffness above e* versus elas- 
tic stiffness, finds expression, we find: 



- p* = 0.640 l(£i - e") , 



p - -o** = 0.624 E«( € i- e **) 

(6.7) 

at u = 0.3, according to Timo shenko ' s calculation. 
It is seen that, according to these formulas, the load 
capacity is rated too favorable, as it should he; because 
Timoshenko uses the energy term without omission while em- 
ploying a Eitz formula, which does not express the actual- 
ly occurring conditions as well as the one used in the 
present report. 

By an extension of (6.1) the cited authors seek to 
account for the fact that the buckles must become deformed 
if the buckling load has been exceeded to a comparatively 
appreciable extent when, as a matter of fact the bulge form 
(6,1) is derived from the energy balance (between bending 
and strain energy) at the instant of buckling, and it is 
evident that this balance must change as the compressive 
forces move toward the sides. The first three authors 
quoted are unanimous in assuming that the profile form of 
the bulge must remain unchanged (w = cos ^ v|>(y)) the 

compressive (x) direction, but that transversely to it, a 
flattening takes place in the middle, and a corresponding 
buckling near the edges. According to that assumption, 
the longitudinal fibers subjected to higher normal compres- 
sion are curved much more, and therefore, are better able, 
to avoid the compression, and the increase in bending energy 
is less than the thus resulting decrease in strain energy, 
Oox particularly assumes the profile form to be built up 
from two sine arcs at the edges and a straight piece in 
the middle, the length of which follows from the minimum 
requirement. Without attacking the other assumptions by 
Cox (ay = T = 0), Yamamoto and Kondo turn against the ar- 
bitrariness of this assumption regarding the buckling form, 
which, in fact, due to the jump in the curvature, is not 
at all satisfactory. And that is the reason for their "ex- 
act" computation of the form of buckle by means of a non- 
linear differential equation. As remarkable as their en- 
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suing mathematical problem undoubtedly is, the labor in- 
■volved is not justified in this case b o cau.s.e : f i rs t , the 
solution of the problem stands and falls with the premises 

CTy (x,y) = T (x,y) = 0; that is, its importance from the 

practical point of view is subordinate; second, this method 
is not worth while even in this particular case because, 
with the help of the continuous formula in the curvatures 
(see report cited at the beginning): 

ttx / Try 3Try\ 

w = cos ^- f f x cos - f 3 cos J (6.8) 

The two parameters f x and f 3 give with the energy 
method the same buckle pattern and the same stress-strain 
law with slide-rule accuracy as Yamamoto and Kondo achieved 
with their exact method. 

Schnadel also employed this same formula (6.8). His 
results differ very little from those obtained without the 
assumption r => 0. The difference in. buckle form (6.8), 
and so from the stress-strain law (5.5) is, moreover, 
slight in the case of. a y ^ 0 if only one deformation of 

the profile x = constant is taken into consideration. 
(This remark was also made by Timoshenko.) 

The relatively good agreement of Cox's formula for 
the apparent width under 100 times exceeded buckling ■ 
crushing, is merely accidental. Because the load capacity 
is considerably underestimated at the very beginning on 
account of the disregarded internal transverse and shear- 
ing stresses, and that necessarily is righted again under 
sufficiently exceeded load, because the buckling form, 
consistently diverging more and more from the actual form, 
mistakingly presents a too high load capacity. 

VI. EXTENSION OE FORMULA (6.1) 

The Apparent Width Under Considerably 
Exceeded Buckling Load 

The aforementioned displacement of the pressure dis- 
tribution under considerably exceeded buckling load (as 
often encountered in practical airplane design) calls for 
a much greater divergence from the buckle form (6.1) than 
afforded in (6.8). For it is a fact that under ten times 
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exceeded "buckling load, for example, the two strips adja- 
cent to the edges carry nearly the full load, while the 
middl-e of the plate takes up little more than the "buck- 
ling stress itself. The conception suggested "by the ap- 
parent width, that the plate falls approximately into 
three zones: the two edge strips of width bm (the appar- 
ent width), and the unloaded middle part are indicative 
of the kind of strain to "be expected. Intermediate buck- 
les must form near the edges since the edge strips, ex- 
actly like the whole plate (at start of phenomenon), have 
the tendency to split more into square panels rather than 
into the very long rectangle*, so that we may put 

_ ttx iry ■ 3ttx / Try _ 3ny\ 

w = fi cos — cos •— - f s cos --— (cos ~ - T| cos ---) 

b b 3 b V b b r . • 

(7.1) 

with the parameters f i , % , and T|, For T] = 1 , for 

example, this formula would have the centroidal axis y = 
0 retain its cosine form unchanged, while toward the 
edge, the middle of each panel (x = 0) develops appreci- 
able counterbuckles under sufficiently great f 3 . With 
proper choice of T| < 1 the exact location of these 
buckles can be more accurately determined and a certain 
amount of flattening obtained even in the middle. 

The calculation itself may be made "by either one of 
two methods: the "exact" method, or the approximate meth- 
od evolved on the formulas (6;3) and (6.4). 

The "exact" method being, as is quickly proved, ex- 
tremely tedious, we give hereinafter, only the approxi- 
mate method and then compare the results. 

It is : 

Xx Xy . _ Xx/ Xx m 3Xy\ A 2tt\ 

w=fi cos -- cos-- - f 3 cos 3 --{cos --- T| cos — -- U = ~-^J 

I) w x 3 = i f x 2 (1-cos X x)(l+cos X y) - I f 3 3 (1-cos 3X x) 
[1 + cos X y - 2T](cos X y+cos 2X y)+Tl 2 (l+cos 3X y)], 



The phenomenon of intermediate buckles above ten times 
exceeded "buckling load, was experimentally observed by 
Lahde (Luf tf ahrtf or schung , vol. 13, 1936, p. 214). 
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•| fifsCcos X z-oo.s 2X x) [(1+oos > x)-Tl(cos X y+cos 2X y)], 



8 

•r ) w. 



a - i f x 2 (1+cos X x) (1-cos X y) + ^ f 3 (l+cos 3X x) 



4 



[1 - cos X y - 6Tl(cos X y-cos 2X y) + 9T| 2 (1-cos 3X y)] 



cos X y+cos 2X x) [(1-cos X y)-3Tl(cos X y-cos 2X )], 



C = 



_ 1 



2b 



w x 2 dx 



= {f, 2 (1 + cos X y) 

8b 2 1 



u 

+ 9f 3 a [(1 + cos X y) -2T) (cos X y + cos 2X y) 

+ Tl a (1 + cos 3X y)]| 

b 

- 1 - |" Wy 2 dy = J* {f, 2 (1 + cos X x) 

o 

+. f 3 a [(1 +.cos 3X x) (1 + 9T1 2 )] 
- 2f x f 3 (cos X x + cos 2X x)| 

{5 v 

b <-2 

+ 81 f 3 4 : | (1 + Tl 4 ) + 6T| 3 - 2TlJ 
18 fi 3 f 3 3 1 "+ Tl a + | (1 - 27)) J } 



1 / C a dy = — n=, J 5 - 4 
b J 64 b <-2 

6 



1 
b 



TT 



D'* dx = -~ * , 

64 b V.2 



f 3 4 . 3 „ 4 . _2 ,s 

12 1 



+ | f 3 4 (i + 9Ti a r 



at : 



0 1 
+ f x 3 f 3 2 [4+2(1 + 9T) 2 )^-.2f 1 3 f 3 | 

Then (6.4) gives the strain portion A of the PE 
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A = G s b 3 (c! 2 - £l -JL. [ fl 2 +9 (1 + T\ 3 ) f 3 3 ] 

C 4 b 



+ C 2 3 + C 2 ^-3 [fi 8 + (1 + 9T1 3 ) f 3 2 ] 



+ -f^-* ^f^^f! 3 f 3 +6f x 2 f 3 2 (^-3T)+6T} 2 ) 

+ f, 4 (123 ~ 162T1 + 513T) 2 + 243 T1 4 ) 1 (?.2 3 ) 

J 

and for the "bending portion X, according to (2.7) 
A = ^~|- [4f x a + f 3 2 (100 + 324 T| s )] 

= G_.iL 3 . [ fi 3 + f;j 2 (25 + 8 i ^2)] (7.2 4 ) 
3 

In view of the high degree of the ensuing equations, 
it is not advisable to attempt the determination of f x , 
f 3 , and T) in relation to ei by means of the three equa- 
tions : 



SA_ = o 
9f x * 



BA_ 
3f 3 



0, 



cA 



(7.3x) 



and one of the conditions of section V for e a by the 
usual method, "but rather to abandon the condition ^ = 0 

for T] , the most unessential of the parameters, and to 

arrive by trial and error at a value of T| yielding the 

lowest possible load capacity. The differentiation, ac- 
cording to f x and f 3 gives: 



- € X - € a + 



TT 



32 b 



2- [12 fi 2 + 6 f 3 2 (11-6T1+12T1 2 ) 



- 9 (1 + n 3 ) e x - (1+9T1 2 ) e 2 



- 6f x f 3 ] + e* = 0 
fi 3 



+ 3^™z ^6f x 2 (11-6T1+12T1 2 ) -2 -j- 

+ 4f 3 2 (123 - 162T1 + 513T) 2 + 243T] 4 ) 

+ (25 + 81T] 2 ) €*=' 0 



(7.3 S ) 
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It also follows from p 3 = 0 (the only case considered 
here), according to (6.3) and (7.2i): 



£s = 32 



^ j^ 4fl 2 + 4 (1 + 971 2 ) f 3 2 ] (7.3 3 ) 



that is: 



y (?.3 4 ) 



€l = € * + —H-s [ 4fl 2 _ 3fl f 3 + f3 a (31-1871+18T| 3 )3 

16 D 

9 (1 + Tl 3 ) €i = (25+81T1 2 ) e* 

+ TT-?TT [fi 2 (31-18TI+18T1 2 ) - f x 3 /f 3 

16 D 

+ f 3 2 ( 244 - 3 24T| + 990T1 2 + 324T] 4 )] 



The elimination of the parameter f x from "both equations 
leave s : 



25 + 81 H 

£ i p 

9 + 9 Tl 2 



= 9 (i + n 2 ) 



4 - 3£ + (31 - 18T1 + 18T1 2 ) V 



(31-18T1+18T| 2 ) - l/t + (244 - 324T} + 990T1 2 + 324T1 4 ) £ 2 

with which tj/e* is readily computed for every T| as 
function of £ = ~- , and inversely £ (j*) ' 

Furthermore, p x is given according to (6.3) and 
(7.2) through the equation 



b/2 

i J w « 

-h/2 



CT, dy = e,. - -4 Cf! 2 +9 (1+Tl a ) f 3 2 ] (7.4 2 ) 

8h 



Inserting Cj. acoording to (7.34) gives: 
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Pl = p* + E f^z l2f x a - 3f x f 3 + (13-18T1) f 3 2 ] (7.4 S ) 
that is, 

r " " r '* S 4 - 6t + (26 - 36Tp f (7.4 4 ) 



Pi 



€ X - e* 2 4 - 3£ + (31 - 18T1 + 18T] 2 ) £ B 

This formula shows conclusively that the apparent stiffness 
decreases at £ > 0 (i.e., "by increasing strain of the 
buckles). Analyzing Vi( € i) as follows, from (7.4 X ) and 
(7.4 4 ) for different values of T\ , it is seen that with 
T) = a the most logical course of the stress-strain curve 
is obtained. Then (7.4 4 ) and (7.4J read: 



Pi - P* 



B 4 - 6£ + 8£ 2 

" 2 4 _ ~3T~ri>67ir~P = 



CP U) 



(7.5) 



€j_ - 4.02 e* 



4 - Zt + 26.5 i s , , . 

= 11.25 ^— ~ -rg- 01 = i) 

26.5 - l/t + 350 I 



(This equation (7.5) represents the sought-for stress- 
strain law.) 

Based upon the exact method, i.e., with consideration 
of the part of the FE originating with the shear stresses, 
we obtain (likewise for T) = ■§) in place of (7.5) the re- 
lations? 9 



Ik 



— -Q * 

_T7* 



E 4 - 6£ + 18.6 £ 3 E 

i i"~ "iI~~3lT8~I^ " 2 9 ( ^ 



€ x - 4.02 e* 



= 11.25 



4 - 3£ + 31.8 V 



> (7.6) 



31.8 - l/t + 350 



It will be observed that the majority of the terms in 
these equations (7.6) are in agreement. 10 The discrepan- 



y The difference in both theories touches only the factor 
of the (mixed) term f j 2 f 3 2 in (7.2 3 ) for, the FE, since 

the shear stresses do not contribute to the other terms 
(a useful check for the "exact" method). 

10 The transverse contraction precisely cancels out in this 
considered case p 2 = 0. 
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cies concern only terms which make themselves felt when 
£ = 0 (that is, for e± ^>€*) is con si derahly departed 
from. They act in the same direction in "both equations: 
there is for equal £ value, according to (7.6) a greater 
€i than according to (7.5) and a greater reduction fac-? 
tor cp ( £) . As it should he, the discrepancy from (5,5) is 

less with the exact method than with the approximation law 
(6.4). 

In figure 2, which gives all three curves, the curves 
for the e-L values diverge not only from their tangents 
hut also from each other more and more. The result of the 
approximation, as expedient as fee ahhreviated calculation 
may he for a first approach, must he received with caution; 
its "better agreement with the Lahde-Wagner test points un- 
der con si derahly exceeded Duckling load (fig. 3) is, of 
course, accidental. Cox's theoretical and experimental 
points are also shown in figure 2, and show the good agree- 
ment of his experimental points with our own results. 

Although the equaiions (7.6.) contain the essentials 
regarding the "behavior of the plate after "buckling, it may 
"be of general interest to give the exact formula for the 
apparent width. Expressing the apparent width with 

V PL = Px * ( 7 * 7 i 

gives with p^ = E e x (= stiffener stress) as extension 
of (5.8): 

i..£ + !«p(t> (l |1) (7.7. 

One may attempt to approximate this curve hy a sim- 
ple analytical expression. In support of von Karman 1 s 
approximate formula (reference 6), it suggests the use of 
a formula of the form of 

In fact, a satisfactory approximation is ohtained with 
A = 0.81 and B = 0.19 within the range of e* j£ €]_ <: 

60 e* or, in other words, with the formula 

h m = 1.54 /i- s + 0.19 T> (7.8) 

J PL 
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Within a small intermediate range the values obtained by 
the formula are too small. 

Another " empirical" approximat e formula reads: 



This formula, while giving values which are a little too 
low for d > 20 e* relative to (7.6) reproduces, how- 
ever, the typical "behavior (inclusive of e x = €*) very 
well, and its marked simplicity recommends it. 

The Lahde-Wagner test points in figure 3 are merely 
"by way of reference, since they pertain to the case of 
nonshifting sides and perfect fixity. To "become "compar- 
able" they are "converted"; that is, the course of the 
points with £i/e* is taken over and only the critical 
value c^j, (which, naturally, is higher for fixation), is 

identified with our reference point e*. On the assump- 
tion that this simple conversion method is permissible , 
the agreement, particularly with the formula (7.9) is 
vory good. 



Translation by J. Vanicr, 
National Advisory Committee 
for Aeronautics. 




(7.9) 



n inn ■iiiiw inn i n imiiiii i urn mi m mi 
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Figs. 1,2,3 



! 




s 



Figure 1.- Rectangular plate under uniform compression on 
two sides. 
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a. Theoretical curveCeq. 7.6) 

b. " » (eq. 7.5) 
approximate method 

c. Theoretical curve(eq. 5.6) 
(tangent to a and b) 

d. Cox's theoretical curve 
Figure 2.- Stress-strain curves 

according to different 
theories for range ei/e*<75(p 1 = 0) 
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a. Theoretical curve^eq. 7.6) 

b. Approximation for a) 
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= 0.811/ + 0. 
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c. Approximation for a) 
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d. Cox's theoretical curve. 

Figure 3.- The apparent width after 

exceeding the buckling 
load p** E-c* — theoretical and ex- 
perimental. 
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